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Carbon monoxide (CO) that is produced in a number of different mammalian tissues is 
now known to have significant effects on the cardiovascular system. These include: (i) 
vasodilation, (ii) changes in heart rate and strength of contractions, and (iii) modulation 
of autonomic nervous system input to both the primary pacemaker and the working 
myocardium. Excessive CO in the environment is toxic and can initiate or mediate life 
threatening cardiac rhythm disturbances. Recent reports link these ventricular arrhythmias 
to an increase in the slowly inactivating, or "late" component of the Na + current in 
the mammalian heart. The main goal of this paper is to explore the basis of this 
pro-arrhythmic capability of CO by incorporating changes in CO-induced ion channel 
activity with intracellular signaling pathways in the mammalian heart. To do this, a quite 
well-documented mathematical model of the action potential and intracellular calcium 
transient in the human ventricular myocyte has been employed. In silico iterations based 
on this model provide a useful first step in illustrating the cellular electrophysiological 
consequences of CO that have been reported from mammalian heart experiments. 
Specifically, when the Grandi et al. model of the human ventricular action potential is 
utilized, and after the Na + and Ca 2+ currents in a single myocyte are modified based on 
the experimental literature, early after-depolarization (EAD) rhythm disturbances appear, 
and important elements of the underlying causes of these EADs are revealed/illustrated. 
Our modified mathematical model of the human ventricular action potential also provides 
a convenient digital platform for designing future experimental work and relating these 
changes in cellular cardiac electrophysiology to emerging clinical and epidemiological data 
on CO toxicity. 

Keywords: carbon monoxide (CO), late sodium current (Iwa-L). action potential (AP), early after-depolarizations 
(EADs), S-nitrosylation 



INTRODUCTION 

It is now well known that carbon monoxide (CO) has significant 
physiological and pathophysiological effects in a number of dif- 
ferent organ systems in mammals (Wu and Wang, 2005). This 
paper considers CO induced ventricular rhythm disturbances 
in the human heart. Recently, attention has been drawn to the 
fact that CO can induce cardiac arrhythmias. This is based on 
a number of studies in animal models and literature regarding 
emergency hospital admissions due to life-threatening cardiac 
rhythm disturbances triggered by increases in ambient levels of 
CO (Goldstein, 2008; Bell et al, 2009; Dallas et al, 2012). More 
generally, substances somewhat similar to CO, including nitric 
oxide (NO), and hydrogen sulfide (H2S) are now often referred 
to in the neurophysiological literature as gaseous transmitters 
(Leffler et al., 2006). These compounds are receiving increased 
attention as (i) mediators or modulators of significant, some- 
times acute pathophysiological challenges, and as (ii) potential 
targets for therapeutic interventions (Motterlini and Otterbein, 
2010). 



The ability of CO to alter electrophysiological activity in the 
heart and nervous system has led to quite broad-ranging exam- 
inations of the ability of CO to regulate specific ion channels 
(Wilkinson and Kemp, 2011). These effects can occur either as 
direct actions of CO, or (and more commonly) CO effects on 
closely related intracellular signaling pathways, often involving 
nitric oxide signaling, metabolism, and/or downstream targets for 
chemical modifications of individual amino acids in functional 
proteins. In the cardiovascular system, specific attention has been 
drawn to S-nitrosylation of, specific ion channel residues (faffrey 
et al, 2001; Haldar and Stamler, 2013). 

In this paper, our focus will be on the cellular electrophysio- 
logical effects of NO and related S-nitrosylation of specific sites 
on ion channels, ion transporters, or directly related signaling 
molecules (Jaffrey et al., 2001; Tamargo et al., 2010; Haldar and 
Stamler, 2013). This is because many cellular effects initiated 
by increased CO levels are mediated by downstream changes in 
NO and ultimately through targeted nitrosylation of defined sites 
(residues). At present, it is known that NO/S-nitrosylation can: 
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(i) increase the background K + current that generates the resting 
potential in heart (Gomez et al., 2009) (ii) activate ATP sensitive 
K + channels in neurons (Kawano et al., 2009), and (iii) reduce a 
major repolarizing current (K v 1.5) in mammalian atria (Nunez 
et al., 2006). In addition, a similar CO and/or NO induced reac- 
tion mechanism, as S-nitrosylation, can modulate (iv) Na + /K + 
ATPase turn over in the hypoxic heart (Yakushev et al., 2012), and 
(v) alter release of Ca 2+ from the sarcoplasmic reticulum (SR) 
by targeting the ryanodine receptor complex in heart and skele- 
tal membranes (Gonzalez et al., 2010; Wang et al., 2010). There 
is also evidence that NO can increase the slowly inactivating or 
late Na + current in nerve, muscle and heart preparations (Ahern 
et al, 2000; Evans and Bielefeldt, 2000). 

Two recent papers (Abramochkin et al, 2011; Dallas et al., 
2012) have reported CO induced changes in electrical activity 
and contractions, as well as induction of cardiac arrhythmias 
in a rat ventricle model. Evidence that a major ion channel 
mediated effect is an augmentation of slow inactivation of the 
cardiac Na + current is provided. These results also emphasize an 
essential involvement of NO as a second messenger of this CO- 
induced effect. CO-induced increases in this late Na + current, 
iNa-L; can significantly lengthen the action potential (AP) and 
result in abnormal electrical activity (early after-depolarization, 
EADs) characterized by repetitive firing, even after application of 
only one stimulus (Dallas et al., 2012). This arrhythmia can be 
reduced substantially following superfusion of this experimental 
preparation with the compound ranolazine (Dallas et al., 2012). 
Ranolazine, originally developed as a coronary vasodilator, is now 
known to be a potent and quite selective inhibitor of the INa-L in 
a number of different mammalian preparations (Makielski and 
Valdivia, 2006; Zaza et al., 2008), and also in human ventricle 
(Moss et al, 2008). Ranolazine effects are significant in the set- 
tings of stable angina and also in heart failure where, in both cases, 
it is known that lN a -L is increased (c.f. Trenor et al., 2012). 

The main goal of our mathematical simulations is to illustrate, 
using a current mathematical model of the AP in human ven- 
tricle, the consequences of a CO induced: (i) increase in Ijsia-L, 

(ii) NO-S-nitrosylation induced changes in Ic a -L> and (iii) this 
combination. For this purpose, CO-induced changes in I^a-L that 
are similar to those reported in the primary experimental data 
(Dallas et al., 2012) were introduced into two different mod- 
els of the ventricular AP and the consequences were explored 
systematically. 

METHODS 

HUMAN VENTRICULAR MYOCYTE MODELS 

To simulate the electrical activity of human ventricular myocytes, 
two relatively recent AP models were evaluated. The Grandi et al. 
model (Grandi et al., 2010) is perhaps the most detailed math- 
ematical model that is available, when judged in terms of its 
comprehensive ionic current portfolio and inclusion of mathe- 
matical expressions for the Ca 2+ transient and Ca 2+ homeostasis 
in the human ventricle. For the purpose of simulating some 
aspects of heart failure we have previously modified and utilized 
the original Grandi et al. model (Trenor et al, 2012). In this paper, 
the slowly inactivating component of lN a , which is denoted lN a -L> 
was reformulated based on experimental data. In the simulations 
presented here, we have made use of this code: specifically the 



Grandi model was employed after (i) modifying the Na + current 
as in (Trenor et al., 2012), (ii) changing the L-type Ca 2+ current, 
Ica-L, or (iii) both. 

Early in this study, some simulations were carried out using 
the latest human ventricular AP model published by O'Hara et al. 
(2011). This model is based on experimental data taken from 140 
healthy human hearts; it encompasses the formulation of 18 ionic 
currents and carrier-mediated fluxes and a detailed formulation 
of steady-state and transient ion concentrations, including intra- 
cellular Ca 2+ transients. This model reproduces many aspects of 
the electrophysiological behavior of human ventricular myocytes 
with high fidelity, and can simulate some AP alterations due to 
drug effects. However, it appears not to be able to generate/ exhibit 
any EADs in response to the reported, CO induced, changes 
in lNa or lea or their combination. We have not explored the 
reasons for this in detail; however, it is likely that the net cur- 
rent at the level of the plateau of the AP is an area for further 
examination/modification. 

SIMULATION OF CO EFFECTS 

As reported in (Dallas et al, 2012), CO can reduce peak tran- 
sient inward lNa by as much as 50%, shift the inactivation curve 
in hyperpolarizing direction, and significantly increase the late 
sodium current INa-L- Thus, in our simulation CO effects on 
the fast component of lNa were modeled by introducing a 50% 
decrease in I^a maximum conductance and a 5 mV shift in the 
hyperpolarizing direction of the two Na + current inactivation 
relationships and j'oo (see Equation 1). In addition, INa-L was 
increased 2-fold compared to its values in the original Grandi 
et al. model (Trenor et al., 2012). 

^oo_shifted = 7oo_shifted = ~~ ~ +71 55 + 5 7 2 ^ ^ 

M + e " 7.43 ) 

V m is the membrane potential. 

CO is known to elevate NO levels (Dallas et al, 2012), and this 
"second messenger" can have important effects on L-type Ca 2+ 
current, Ic a -L> in the mammalian (ferret) heart, as documented 
in the detailed studies of Campbell et al. (1996). These investi- 
gators reported a significant increase in Ic a -L (30-50%) under 
conditions of NO induced S-nitrosylation of this Ca 2+ channel 
a-subunit. In addition, this paper also reported a nitrosylation 
induced change in Ic a -L gating — specifically a small but signifi- 
cant change in the voltage dependent activation relationship — an 
approximately 6.5 mV shift in the hyperpolarizing direction. To 
reproduce these effects using the Grandi et al. model (Grandi 
et al., 2010), we have made a number of changes in these param- 
eters. The most favorable/realistic results were obtained when the 
maximum conductance for Ic a -L was increased by 20% and its 
activation gate was shifted 3 mV in the hyperpolarizing direction 
(see Equation 2). The changes applied were slightly smaller than 
the ones reported experimentally, but were sufficient to trigger 
EADs, as shown in Figure 2. 

^oc_shifted = ~ y m + 5+ 3 (2) 
1 + e *~ 

V m is the membrane potential. 
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SIMULATION OF RANOLAZINE EFFECTS 

The effects of ranolazine, at an assumed plasma level of 5-10 u,M, 
were simulated by decreasing the maximum iNa-L conductance by 
50% (see Figure 4) (cf. Trenor et al, 2013). This maneuver was 
based on the fact that the recommended adult dosage level for 
ranolazine is in the 3-8 [lM range (Belardinelli et al., 2006). 

STIMULATION PROTOCOLS 

AP simulations were conducted at a stimulation rate of 1 Hz. 
Measurements were taken on stimulated output data only after 
achieving steady-state conditions. 

NUMERICAL IMPLEMENTATION 

All model equations were taken from Grandi et al. (2010), and 
were implemented in Matlab (Mathworks Inc., Natick, MA, 
USA). Differential equations were solved numerically using a 
variable order solver (odel5s) (Shampine and Reichelt, 1997). 
As indicated, some simulations were performed using O'Hara 
et al. model, which was also implemented in Matlab. These model 
equations were downloaded from http://rudylab.wustl.edu, and 
the rapid integration methods provided in the Supplemental 
Materials from O'Hara et al. (2011) and one of our previous 
papers (Trenor et al., 2012) were used. 

RESULTS 

The initial set of computations was done for the purpose of 
illustrating the consequences of the CO induced changes in mam- 
malian heart Na + current, I^a that were described in a primary 
experimental study. This paper (Dallas et al., 2012) reported that 
CO acting through an intracellular NO-mediated signaling cas- 
cade can (i) reduce peak lN a by as much as 50% while also (ii) 
shifting the inactivation curve for the large transient compo- 
nent of 1^ in hyperpolarizing direction, and (iii) significantly 
increasing the slowly inactivating component of 1^ (which is 
denoted In 3 -l)- 

The two superimposed human ventricular APs in Figure 1A 
illustrate the control AP waveform at 1 Hz, and the AP elicited 
following the changes in lN a and In 3 -l described above. For these 
computations, the steady-state inactivation curve was shifted in 
the hyperpolarizing direction by 5mV and INa-L was increased 
2-fold. Figure IB shows the changes in INa-L; Figure 1C illustrates 
the indirect effects on Ica-LJ that is, the AP waveform lengthened, 
and this changed the Ica-L! and Figure ID shows the computed 
intracellular Ca 2+ transient, [Ca 2+ ];. 

In summary the reported changes in lN a and lN a -L> when 
combined, result in predictable decreases in AP rate-of-rise 
(not shown) and small increase in AP duration. However, 
these changes fail to elicit spontaneous firing (EADs on DADs) 
or produce any correlates of arrythmogenesis. Although this 
pattern of results was somewhat unexpected, it is impor- 
tant to recall that the primary experimental data on CO 
were obtained in rat hearts (Ahern et al., 2000; Dallas et al., 
2012), and not from human ventricular tissue or myocytes. 
We therefore, continued to attempt to illustrate and fur- 
ther understand the electrophysiological effects of CO on 
the human ventricular myocardium by considering additional 
variables. 
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FIGURE 1 | Effects of carbon monoxide (CO) induced changes in Na + 
current on the action potential in a human ventricular myocyte. In this 
and all subsequent Figures, the Grandi et al. mathematical model of the 
human ventricular action potential was employed (see Methods). Control 
simulations are shown as black traces and simulations generated by 
selected changes in model parameters are shown in red. Panel (A) consists 
of control action potentials (black, 1 Hz., steady state) together with 
superimposed action potentials (red) that were elicited by the same 
stimulus parameters after: (i) reducing peak Na+ current by 50%, (ii) shifting 
the steady state inactivation relationship 5mV in the hyperpolarizing 
direction, and (iii) altering the mathematical descriptors that govern the 
slowly inactivating or late sodium current, I Na-Lr such that its size was 
increased 2-fold. As expected, these changes reduced excitability and 
decreased the dV/dt of the action potential upstroke as well as lengthening 
the AP at both APD-30 and APD-90. l Na -L records are shown in Panel (B). In 
Panels (C,D) the corresponding L-type Ca 2+ current, l Ca _ L , and intracellular 
calcium transient, [Ca 2+ ]i are illustrated. Note that these changes are 
indirect, i.e., they result from the changes in APD waveform and the 
intrinsic biophysical properties of Ica-L and tne Ca 2+ -induced Ca 2+ release 
(CICR) mechanisms. 
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In fact, studies of NO-mediated effects in heterologous expres- 
sion systems utilizing the alpha subunit of the human heart Na + 
channel, Na v 1.5 (Ueda et al., 2008), reveal a somewhat differ- 
ent pattern of results. These include (i) an (approximately 30%) 
increase in peak inward Na + current, a (ii) small hyperpolarizing 
shift in the steady-state inactivation curve for I^a, (hi) a marked 
enhancement (approximately 2-fold) in iNa-L and (iv) an essential 
NO-induced S-nitrosylation of a defined residue on Na v 1.5. These 
changes, when introduced into the Grandi model, also failed to 
result in increased arrhythmogenesis (results not shown). 

We were aware of the detailed studies of Campbell et al. 
(1996) on the effects of NO on Ic a -L> in the mammalian (fer- 
ret) heart. These investigators reported a significant increase in 
Ica-L (30-50%) under conditions of NO induced S-nitrosylation 
of this Ca 2+ channel a-subunit. Importantly, this paper also 
reported a nitrosylation-induced change in Ic a -L gating, specifi- 
cally, a small but significant shift in the voltage dependent acti- 
vation relationship — approximately 6 mV in the hyperpolarizing 
direction. The superimposed APs in Figure 2 illustrate the con- 
sequences of introducing these changes into the original Grandi 
et al. model (Grandi et al., 2010). Panel (A) shows the AP wave- 
forms (black, baseline, or control) (red, following Ic a -L changes). 
INa-L and Ica-L records are shown in Panels (B,C), respectively. 
The L-type Ca 2+ current is shown in Panel (C). As shown in 
Panel (D), these changes in Ic a -L also resulted in a large increase 
in [Ca 2+ ];, and the [Ca 2+ ]; exhibited pronounced frequency 
dependence or alternans (see Discussion and Clark et al., 1996; 
Bouchard et al., 2004). 

Detailed inspection of Figure 2 and comparison of this pattern 
of results with those in Figure 1 reveal (Panel B) that INa-L now 
shows pronounced changes in peak amplitude, even in response 
to the fixed stimulus rate at 1 Hz. Presumably, this arises from 
the markedly prolonged AP and intrinsic voltage-dependence of 
reactivation of INa-L (see Discussion). 

The third set of computations in this study combined the 
reported CO induced changes in In s and I^a-L (Figure 1), with 
those for Ica-L (Figure 2). The resulting very significant changes 
in AP duration (Panel A), INa-L (Panel B), Ica-L (Panel C), and 
[Ca 2+ ]; (Panel D) are shown in Figure 3. This pattern of results 
was similar to, but not identical with Figure 2. In particular, the 
majority of the APs elicited at a physiological stimulus rate (1 Hz) 
exhibited marked prolongation (approximately 100%), with the 
appearance of an EAD late in the prolonged plateau phase. Note 
that under these starting conditions, the combined inward cur- 
rents due to INa-L and Ica-L "held" the membrane potential of the 
plateau near OmV in contrast with the approximately — lOmV 
plateau level in Figure 2. 

The pattern of changes shown in Figure 3 is not surprising. It 
is now well known that abnormal, slow depolarizations, which are 
denoted EADs, can be elicited by changes in Ic a -L (Marban et al., 
1986; January and Riddle, 1989, for review see Clusin, 2003). In 
fact, in rabbit ventricle a recent report from the Weiss laboratory 
(Madhvani et al., 20 1 1 ) has illustrated in detail ways in which even 
very small changes in the voltage dependence of the activation of 
Ica-L can give rise to EADs. This is mainly, but not entirely due 
to the fact that there is an increased inward current at or very 
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FIGURE 2 | Effects of CO-induced changes in L-type Ca 2+ current on 
the human ventricular myocyte. Based on the work of Campbell et al., 
Ica-L was increased 50% and the voltage dependence for the activation 
curve was shifted in the depolarizing direction by 3mV. As shown in Panel 
(A) these changes resulted in a lengthening of the action potential plateau 
and appearance of associated oscillatory potentials or early after 
depolarizations (EADs). Panel (B) shows the changes in l^a-L- The longer 
APD results in a protracted l Na . L ; the alterations in amplitude of this 
current are due to the voltage dependence of its reactivation. As shown in 
Panel (C), these changes in the biophysical parameters that regulate lc a -L 
resulted in a change in peak amplitude and perhaps more importantly a 
significant alteration in its time course of inactivation. Ic a -L also reactivates 
under these conditions. Panel (D) shows the control Ca 2+ transient and 
the computed changes in [Ca 2+ ]i that occur as a consequence of 
alterations in the lc a -L- 
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FIGURE 3 | Changes in the human ventricular action potential as a 
consequence of CO-induced changes in both lfj a and lc a -L- The layout 
and color-coding of this Figure is identical to that described for 
Figures 1, 2. As shown in Panel (A) when the CO-induced changes in 
Na + and Ca 2+ currents are combined there is marked lengthening of the 
action potential and a significantly increased incidence and frequency of 
EADs occurring during the somewhat depolarized plateau phase of the 
action potential. l Na _ L records are shown in Panel (B). Panels (C,D) show 
the changes in lc a -i_ and [Ca 2+ ]j. Note that the changes in lca-L differ 
from those in Figure 2 as a result of the alterations in action potential 
waveform caused by the combined changes in l Na and Irja-L- 
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FIGURE 4 | Estimation and illustration of the effects of ranolazine on 
CO-induced changes in the action potential and underlying ln a -L- 
'ca-L- a "d Na + /Ca 2+ exchange currents (Incx) human ventricle. The 

effects of ranolazine (shown in black) were modeled as a selective 50% 
reduction in the l Na . L . Note that this reduces action potential, duration of 
every 3rd action potential, as well as decreasing the incidence and the 
frequency of EADs. However, these changes are complex. APD (Panel A) 
as well as |Ca 2+ ]i and Incx records all show variations or alternans, even 
at this fixed stimulus rate (1 Hz). Panel (C) illustrates the corresponding 
ranolazine-induced changes in the Incx- Note that in the action potential 
traces showing a ranolazine-induced shortening of APD, there is no 
spontaneous inward l NCX current (see also Figure 5 and Results). The 
arrows denote the "extra" EADs (Panel A), and underlying [Ca 2+ ]i 
transient (Panel B), and Incx (Panel C) records. 
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near the range of membrane potentials of the plateau of the nor- 
mal or prolonged AP. The underlying biophysical mechanism(s) 
involves the non-linear interactions of the small overlapping Ca 2+ 
channel — and Na + /Ca 2+ — exchanger mediated currents in this 
range of membrane potentials (Madhvani et al., 2011, for review 
see Fink et al., 2006). The altered AP waveform and related pro- 
or anti-arrhythmic effects are sometimes described in terms of 
an altered "repolarization reserve" at defined time-points of early 
or late repolarization (Wu et al., 2009; Varro and Baczko, 2011; 
Trenor et al, 2013). 

If environmental or tissue-derived CO can induce life- 
threatening ventricular rhythm disturbances through a mech- 
anism involving enhancement of iNa-L> it is plausible that, as 
mentioned in the Introduction, ranolazine could be effective in 
Critical Care settings (Dallas et al, 2012, cf. Belardinelli et al., 
2006; Bell et al, 2009). The mathematical simulations shown 
in Figure 4 address this possibility by illustrating the effects of 
ranolazine "in silico" These computations were carried out under 
starting conditions (i.e., an initial parameter set) almost identi- 
cal to those for Figure 3. However, for these computations it was 
assumed that ranolazine actions would be the equivalent of a 50% 
reduction of the CO induced increase in I^a-L (see Figure 1). The 
results show that if ranolazine (at an assumed plasma level of 
5 (xM) did block approximately 50% of lN a -L> it could shorten the 
AP elicited by a 1 Hz train of stimuli. Note, however, that this par- 
tial block of iNa-L reduced APD only (approximately) every 3rd 
AP in this 1 Hz AP train. An additional effect is also apparent: 
reducing lN a -L changes the (i) waveform and (ii) the most negative 
membrane potential in some repolarizing components within this 
train of AP's. This primary "effect" appears to be linked function- 
ally to the approximately normal repolarization of the next AP. 
In summary, when ranolazine is introduced in silico, the result- 
ing AP waveforms show complex alternating patterns. However, 
there is an increased number of (i) near normal AP waveforms 
and (ii) monotonic [Ca 2+ ]; records. The incidence of EADs is also 
reduced, consistent with the previously described actions of this 
drug in the settings of free radical challenge, genetic abnormali- 
ties in Na + channel biophysics, hypoxia, or global ischemia/heart 
failure (Fink et al, 2006; Trenor et al, 2012). 

The simulations in Figure 5 address and attempt to illustrate 
some of the underlying causative factors for EAD generation 
within the framework/limitations of this model of the human 
ventricle AP. EAD generation requires the membrane poten- 
tial to be relatively depolarized, so that L-type Ca 2+ channels 
can be activated and contribute a small but significant inward 
Ca 2+ current. We note that the main interacting (or overlap- 
ping) ion channel mediated currents at the plateau of the AP are: 
(i) lNa-L) (ii) Ica-L> (iH) Iri> a time-independent or background 
inwardly rectifying current, and (iv) the delayed rectifier K + cur- 
rent, Ik-r (Trenor et al., 2013). In the presence of an enhanced 
lNa-L the AP plateau is somewhat depolarized and the AP tends 
to be prolonged. In addition, under these conditions Incx gener- 
ated by the Na + /Ca 2+ exchange mechanism can be functionally 
important. A key question is: are the pronounced changes in 
[Ca 2 +]i the cause or the consequence of the EAD/CO-induced 
pro-arrhythmic substrate. A plausible answer can be obtained by 
consideration of the results in Panels (A-F) of Figure 5. 



Panel (A) shows two superimposed pairs of APs, each elicited 
at steady state (1 Hz) before (black) and after (red) introducing 
the changes shown in Figure 3. Panels (B,C) show the corre- 
sponding INa-L and Ica-L records. Although this information 
is redundant with that in Figure 3: here, it provides a means 
for direct comparison with the new data: Panel (D), [Ca 2+ ]; 
records; Panel (E), Incx; and Panel (F), the Ca 2+ content of 
the SR. 

Taken together, these results suggest that in this model the 
CO/NO-induced changes in Ic a -L are critical. With each stim- 
ulus, a large "extra" Ca 2+ influx prolongs the plateau and 
results in enhanced Ca 2+ release and a much larger increase in 
[Ca 2+ ];. This change in [Ca 2+ ]j can stimulate the Incx- Perhaps 
more importantly, however, this time-dependent and cumula- 
tive change in myocyte Ca 2+ homeostasis results in markedly 
increased Ca 2+ in the SR. Eventually, this causes anomalous Ca 2+ 
release from the SR, and the Na + /Ca 2+ exchanger (denoted by 
arrows) responds by generating an inward current that gives rise 
to the EAD. Under these conditions, approximately every sec- 
ond AP is associated with spontaneous release at the steady-state 
stimulus rate of 1 Hz. 

This scheme/scenario seems plausible, since similar patterns of 
changes have been described experimentally. However, it may not 
be readily apparent how ranolazine, acting only on lNa-L could 
be anti- arrhythmic (see Discussion). In brief, the answer appears 
to be that ranolazine, by blocking lNa-L, can shorten the AP. This 
reduction in APD results in reduced Ca 2+ influx through Ic a -L 
during each AP and thus, there is less Ca 2+ -loading of the SR. 
The end result is that SR Ca 2+ content does not reach the abnor- 
mal levels required for anomalous and spontaneous Ca 2+ release. 
Hence, there is no intracellular Ca 2+ "trigger" for the electro- 
genic Incx> that has an important role in initiating the EAD in 
this model (Houser, 2000; Clusin, 2003). 

DISCUSSION 

SUMMARY OF MAIN FINDINGS 

Our mathematical modeling illustrates one set of conditions 
under which human ventricle myocytes can respond to a sig- 
nificant CO induced challenge (environmental or intrinsic) by 
generating an aberrant cellular electrophysiological pattern of 
responses, denoted EADs. 

When using a modified Grandi et al. modeling framework 
(Grandi et al., 2010), the reported (Dallas et al., 2012) changes 
in Na + current following CO exposure were not sufficient to pro- 
vide the required pro-arrhythmic substrate. Additional changes 
in lNa-L parameters may have given rise to EAD waveforms. 
However, we evaluated the possibility that other ion channel 
targets were involved in EAD initiation. Accordingly, our analy- 
ses of CO-induced EADs include changes in lN a -L and Ica-L- As 
shown in Figure 2, even very small changes in Ica-L can produce 
EAD generation. Although this pattern of results has not been 
reported when CO is used as the "toxic" stimulus; it is known 
that NO is a primary mediator of CO-induced pro-arrhythmic 
events in the heart (Dallas et al, 2012). Campbell et al. (1996) 
have established that NO driven S-nitrosylation can modify Ica-L- 
We adopted this approach and were guided by this comprehensive 
biophysical/molecular pharmacological data set. 
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FIGURE 5 | Diagram illustrating a plausible mechanism for generation 
of EADs in our modified Grandi model of the human ventricular 
action potential. Panels (A through C) each consist of a pair of 
superimposed records for the AP l^a-L. and lc a -L records, respectively. 
Those in black correspond to control or baseline conditions, and the 
CO-induced changes are shown in red. Important new information is 



illustrated in the corresponding (i) l N cx currents, Panel (D); (ii) [Ca 2+ ]| data, 
Panel (E); and the current generated by Ca 2+ release from the 
sarcoplasmic reticulum l re |, Panel (F). The arrows denote the EAD 
(Panel A) and underlying spontaneous Ca 2+ release from the SR (Panel F), 
and Incx current driven by the increase in [Ca 2+ ]j (Panel E). Refer to 
Results and Discussion sections for further details and explanation. 



It is interesting that in spite of an apparent "double hit" being 
required for the initiation of EADs, in silico, reduction of only one 
of these contributing factors, a 50% decrease of lNa-L> resulted in a 
reduction in the EAD incidence (see Figure 4). This is consistent 
with the experimental observation that selective pharmacological 
inhibition of I^a-L can serve as an antidote to CO-induced rhythm 
disturbances in rat ventricle (Dallas et al, 2012). 

MECHANISTIC INSIGHTS 

Some aspects of the cellular/biophysical mechanism(s) respon- 
sible for EAD generation can be deduced/confirmed by careful 
inspection of the new data from our simulations. EADs are 
slow, small depolarizations that arise in the ventricle (and in 
Purkinje tissue) within a transmembrane voltage range corre- 
sponding to the plateau of the ventricular AP (approximately — 10 



to +10 mV). A net inward current is required, and the work 
from the Marban (Marban et al., 1986) and January (January 
and Riddle, 1989) Laboratories has established that there is an 
obligatory requirement for activation of Ica-L- In fact, Ic a -L> in 
most cases, must reactivate . That is, a small fraction of the same 
Ica-L channels that have opened previously in the AP to pro- 
duce the distinctive plateau of the human ventricular AP open 
again. This produces a "pedestal" or apparently non-inactivating 
Ica-L component (Corrias et al., 2011). Our simulations are 
interesting in this regard. Detailed inspection of Figures 2-4 
in fact show that EAD generation in this model occurs at a 
time when the Ic a -L pedestal shows a small, transient decrease. 
In the Grandi et al. model, this is due to the specific math- 
ematical expression of [Ca 2+ ]; dependent inactivation that is 
employed. 
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However, the non-inactivating component of Ica-L and the 
inward current due to iNa-L markedly prolong the AP, and this 
change in APD waveform results in a significant increase in Ca 2+ 
influx since Ica-L does not fully inactivate. This extra Ca 2+ influx, 
results in a time-dependent and marked increase in SR Ca 2+ 
content. Eventually this causes a spontaneous and intermittent 
(anomalous) pattern of SR Ca 2+ release. Each of these increases in 
[Ca 2+ ]; gives rise to a significant "new" inward current generated 
by the sarcolemmal Na + /Ca 2+ exchanger mechanism. 

It is worth noting that during CO challenge (intoxication and 
recovery) or in any setting involving stress, the known effects of 
stimulation of the sympathetic nervous system would be expected 
to augment both Na + and Ca 2+ currents. These changes, alone or 
in combination, would enhance Ic a -L and perhaps iNa-L> both of 
which would contribute a net inward current at the plateau level 
of the ventricular AR The Koren Laboratory (Liu et al., 2012) has 
recently developed a robust model of EAD induced arrhythmias 
based on genetically targeted malfunction of distinct repolarizing 
K + currents in adult rabbit ventricle. Their findings provide inter- 
esting insights into underlying mechanisms of pro -arrhythmic 
impulse generation at relatively depolarized (AP plateau) levels 
in the ventricular myocardium. Specifically, they have reported 
enhanced EAD activity either (i) following augmentation of Ic a -L 
with isoproterenol, or (ii) a reduction of repolarizing K + cur- 
rents) during hypokalemic conditions. Both of these maneuvers 
can alter the net current at the plateau of the AP and thus, 
change the repolarization reserve during the early repolariza- 
tion phase of the AP (cf. Fink et al, 2006; Sarkar and Sobie, 
2011; Trenor et al, 2012). Further insights into the ionic mech- 
anisms responsible for EAD initiation and maintenance requires 
definition of, and agreement on, the range of membrane poten- 
tials at which those small oscillating events occur (Gaur et al., 
2009). 

FUTURE DIRECTIONS 

As indicated in the Introduction, there are a number of additional 
targets or sites, at which CO or CO-generated signaling molecules 
are known to act. Perhaps the most important of these are the K + 
currents that are expressed throughout the mammalian heart, and 
have important roles in setting the resting potential and modu- 
lating repolarization. The Tamargo Group (Gomez et al., 2009) 
has shown that the K + conductance which underlies the rest- 
ing potential in atrium and ventricle, K; r 2.1 or Iki, is a target 
for S-nitrosylation. This chemical modification of this current 
results in increased Iki- Although Iki is non-linear, it increases 
the repolarization reserve during final repolarization, and may 
cause a small (1-3 mV) hyperpolarization of the resting potential 
(cf. Fink et al, 2006). More detailed studies of CO effects should 
include this target of CO-mediated channel effects. 

Tamargo et al. (Nunez et al., 2006) have also reported that 
a prominent repolarizing current in the atrium of most mam- 
malian species, K v 1.5, can be altered by S-nitrosylation. K v 1.5 is 
decreased by S-nitrosylation, and this would be pro-arrhythmic. 
We note however, that K v 1.5 is a small current in human ventricle 
(cf. Fink et al, 2006; Grandi et al, 2010). 

As noted in the Introduction, the Makielski Group (Ueda et al., 
2008) has reported that NO and an associated S-nitrosylation 



modification of lN a can result in a significantly enhanced iNa-L- 
In fact, however, their investigation focused on NO effects on lN a 
in the setting of an identified mutation in syntrophin that results 
in one variant of a pattern of ventricular rhythm disturbances, 
denoted the long Q-T syndrome (cf. Moss et al., 2008; Thomas 
et al, 2008). 

It has been proposed (Dallas et al, 2012) that considera- 
tion should be given to utilization of ranolazine or ranolazine- 
like compounds in clinical settings (Emergency or Critical 
Care Medicine) that involve documented CO poisoning and 
involve ventricular rhythm disturbances. Improved clinical inter- 
ventions directed toward minimizing cardiac risk during the 
required recovery period from CO intoxication are much needed 
(Goldstein, 2008; Motterlini and Otterbein, 2010). In prin- 
ciple, and as illustrated, ranolazine or ranolazine-like com- 
pounds could be considered for this purpose (Belardinelli 
et al, 2006; Kinobe et al, 2008; Saint, 2008; Wu et al, 2009). 
Ranolazine quite selectively targets iNa-L- However, the net cur- 
rent change that is the direct result of ranolazine interaction 
with Na + channel, is integrated in the overall cardiac elec- 
trophysiological system (biophysical mechanism) for mainte- 
nance and generation of the plateau of the AP, in a way that 
involves complex, often non-linear interactions of a number 
of different currents (Fink et al., 2006). Thus, a monotonic 
dose-response relationship for ranolazine effects is not antic- 
ipated in this type of clinical/therapeutic setting/application 
and individual differences in the electrophysiological sub- 
strate are expected to/need to continue to be considered 
(Sarkar etal, 2012). 

The computations presented in this paper have all been carried 
out assuming that the CO target is a single ventricular myocyte. 
Although this is a valid starting point, it is known that the isolated 
single human ventricular myocyte has electrophysiological char- 
acteristics somewhat different than the same myocyte functionally 
joined to other myocytes through gap junction or "placed" in a 
syncytium consisting of many other similar or identical myocytes. 
Often, the AP waveform in the syncytium is slightly but sig- 
nificantly different than that in a single cell. For example, the 
plateau "height" is less depolarized in the syncytium. Intracellular 
current flow, and other factors, such as intrinsic electrophysio- 
logical heterogeneity of the myocytes, can also result in the AP 
of a myocyte within a syncytium having a more stable resting 
potential together with much less intrinsic action APD variabil- 
ity (Huelsing et al., 2000; Zaniboni et al, 2000; Spitzer et al., 
2006). In addition, it is now known that in both the ventricu- 
lar and the atrial myocardium the fibroblast/myofibroblast cell 
population (Wu and Wang, 2005) may produce significant elec- 
trotonic effects (Maleckar et al, 2009, cf. Baudino et al, 2006). 
It is possible that either CO or its immediate downstream signal- 
ing molecules could alter fibroblast electrophysiology/paracrine 
function. One or more of these factors may alter the pattern of 
response which we have observed in our in silico investigation of 
CO effects. 

LIMITATIONS 

We are aware that as presented our analysis of CO effects on the 
human ventricle is incomplete. Primary reasons for this include: 
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1. The computational platform that we have used is comprehen- 
sive and has been tested in terms of its ability to illustrate 
the functional roles of many of the individual ion channel 
currents that underlie the AP. However, it and most other 
such models have limitations. In this paper, we illustrate 
pronounced effects on AP waveforms and pro-arrhythmic pat- 
terns of response that arise from very small current changes 
that would be consistent with CO effects. The Grandi model 
(Grandi et al, 2010) can be relied on as a computational 
platform for human ventricular myocytes when exploring 
problems that can be addressed with semi-quantitative end- 
points. However, caution is warranted when the input data 
suggests only e.g., a 1-3 mV shift in an activation parameter 
(cf. Terkildsen et al., 2008). 

2. It is apparent that simulation of CO-mediated electrophysio- 
logical effects requires detailed understanding of each major 
component of Ca 2+ homeostasis in the human ventricular 
myocyte. Recent publications show that S-nitrosylation can 
alter ryanodine receptor function and change Ca 2+ -induced 
Ca 2 + release (Ahern et al, 2000; Wang et al, 2010; Cutler 
et al, 2012). Our work does not address this possibility. 
In fact, in the Grandi et al. modeling environment (Grandi 
et al., 2010) detailed studies of the consequences of changes 
in Ca 2+ -induced Ca 2+ release (CICR) are not possible since 
many of the known subcellular details of CICR release are 
not included in this computational platform. After this model 
is developed further and e.g., the existing electrophysiologi- 
cal components are integrated with known properties of E-C 
coupling and mechanical activation (Tao et al., 2011), further 
insights into CO-mediated effects may be able to be obtained 
and Ca 2+ -pump (SERCA 2) mediated relaxation. 

3. Detailed studies of CO-mediated alterations of mammalian 
cardiac electrophysiological phenomena are relatively recent. 
At present, there is incomplete information on clinically rele- 
vant plasma levels, interspecies variations, and specific intra- 
cellular signaling cascades. This information is essential for 
more extensive simulations that could address cellular or 
molecular mechanisms for CO-lN a -NO interactions (Ueda 
et al., 2008), including those arising from ion channel muta- 
tions (Thomas et al., 2008) or changes in intercellular cou- 
pling (Haas and Landisman, 2011), identification or selection 
of therapeutic agents, and improvement of Systems Biology 
(Winslow et al., 2011) approaches to human cardiac sudden 
death in humans. 

4. Our work is not intended to detract from or replace well- 
established principles of CO generation or action (Maines, 
1997) based on the heme oxygenase system. However, two 
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